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Abstract 



Abstract — The acoustic coupling between two air bubbles immersed in water 

is clearly demonstrated. The system is acoustically forced, and its response is 

detected. The experimental results confirm that both theoretically predicted 

eigenmodes, respectively symmetrical and antisymmetrical, do exist. Their 

frequencies, measured as a function of the bubbles spacing, follow theoretical 

estimations within a 10% accuracy. 
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1. INTRODUCTION 



Bubbles play an important role in the sound propagation in everyday life liquids. For 
example, the murmur of the brooks essentially originates, as first suggested by Bragg Jl]||, in 
the oscillations of air bubbles captured and dragged along by the water. The so-called "hot 
chocolate effect", namely the rising of sound pitch when one repeatedly taps the bottom of 
the mug in which some instant coffee or chocolate is being dissolved, is explained by the 
releasing into the water of the tiny bubbles trapped in the powder @,|j]. Bubble dynamics 
and acoustic properties of liquids containing a large number of bubbles have been widely 
studied for a long time p|-|lO|]. A very good and complete review in this domain was achieved 



by Leighton [11]. Inter alia, the problem of the interaction of two neighbouring bubbles has 



been discussed using fluid dynamics tools [12] or the acoustic-electrostatic analogy [13 



Moreover, the free oscillations of a system of two (and even three) air cavities formed in a 
metal plate lying on a water surface have been theoretically and experimentally investigated 
in detail (including cubic nonlinearities) ||14|| . The aim of the present article is to present a 
simple, readily reproducible, experimental study of the forced oscillation regime of a two-air 
bubble system in water. We begin with a short introductory theory in which we show that 
the two-bubble system is mostly equivalent to a set of two magnetically coupled electric 
circuits. 

2. THEORETICAL MODEL 

An air bubble in water will be considered as a perfect sphereQof radius R(t) = Ro + £(t), 
with variation £ much smaller than equilibrium value Rq. It can be shown that £(t) oscillates 



1 The correction to the Minnaert angular frequency due to deviation from the spherical shape can 
be shown to be negligible Jl5|— 117|] . 



with Minnaert's angular frequency uj = y37P /p Po) where 7 is the specific heat ratio 
Cp/C v of air, and Po and po respectively stand for the equilibrium pressure?] and mass 
density of water. This oscillation is damped through several mechanisms: of course the 
acoustic radiation damping (thanks to which the bubble noise is audible), but also the 
viscous and thermal dampings [H,0. We will neglect, in the following simplified theory, the 
last two ones. Moreover, allowing for the typical IKHz acoustic frequency and 1mm bubble 
size we deal with in our experiment, we will neglect any sound propagation in the enclosed 
air. We thus deliberately restrict the present study to the (radial) fundamental resonance 
of the air bubble-water system. 



2.1 One-bubble free oscillation 

Let us consider one bubble with radius Rq immersed in an infinite volume of water at 
equilibrium pressure Po- Let P(f, t) be the actual pressure at site r and instant t. The extra 
pressure p(r, t) is defined as P(f, t) — P . According to Minnaert's assumption, the enclosed 
air undergoes isentropic transformations and its (extra) pressure p(t) is homogeneous inside 
the bubble. Then, neglecting air's inertia as well as the air-water surface tension, p(t) and 
the radius variation are linked by: 

On the other hand, it can be easily shown that (extra) pressure p(r, t) at distance r from 
the center of the bubble follows a d'Alembert-like ID equation, the solution of which exactly 
reads, for r > R : 



2 The pressure difference accross the bubble boundary due to air-water surface tension is about 
1%Pq for a typical radius of 1mm (see [|| eq. (65b)) and will be neglected: Po is also the equilibrium 
pressure of enclosed air. 



p(r, t) = -p Rl 



Rn 



Rn 
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where c is the sound velocity in water. If the acoustic wavelength A is much larger than 
r (i.e., under the circumstances, if condition ruin/c <C 1 is fulfilled), then p(r,t) can be 
approximated by: 



P(r, t) ~ p R 



Po-Rq 



r c 



r c 

Then, equalling p(t) in eq. (|l|) to p(i?o, i) in eq. (H), one gets, all calculations carried out: 



(3) 



,# + u%Rq ^ , 37P0 



^ + 



(4) 



c po^o 

which is the well-known differential equation of a weaklyf] damped ID harmonic oscillator. 



2.2 Two-bubble free oscillation 



Let us now add a second bubble, with the same (equilibrium) radius Rn, at a distance 
d apart from the first one. Let fl (i — 1,2) be the (equilibrium) position of the i th bubble 
center, £j(t) its radius variation, pi(t) the (inner) extra pressure of the enclosed air, and 
Pi(r,t) (resp. Ui(r,t) the would-be (outer)extra pressure (resp. displacement with respect 
to equilibrium) at point r and instant t in the water medium if bubble % was alone. Then, 
allowing for the superposition principle for small displacements, we assume that overall water 
extra pressure and displacement respectively read: 



p(r,t) =p 1 (r,t) +p 2 (r,t) 
u(r, t) = Mi (f, t) + u 2 (r, t) 



(5) 
(6) 



3 Ratio T ra d/ujQ = LOnRn/c is actually assumed to be much smaller than unity, as a consequence 
of the the validity condition of eq. (g) . 
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with, of course, Pi(t) and £j(t) still linked by eq.(|l|). On the surface of the first bubble: 
ri = \r — fi\ = Rq, r 2 = \r — f 2 \ ~ <i, and p(r, t) = pi(t). A similar constraint is required on 
the surface of the second bubble, where r\ ~ d and r 2 = Rq. If the bubble spacing d is much 
smaller than A (i.e. u Q d/c <C 1)0, then eq. (||) is available and we finally get the following 
pair of coupled motion equations: 

el' + a& + r md (e; + q + ^ 2 6 = o (7) 
ad' + e 2 ' + r md (ei + ei) + ^06 = (8) 

where a = Ro/d (< 0.5) can be regarded as a dimensionless coupling constant. Observe, 
by the way, that if double condition: Rq -C d A is fulfilled, eqs. (0) and (H) are available 
(with a ~ 0), and dynamic variables £j are still coupled by radiation damping, since the 
dissipation terms do not involve a. 

Defining symmetrical and antisymmetrical normal variables <p s (t) and <f) a {t) as respec- 
tively the sum and the difference of £i(t) and £2^), we get the uncoupled equations system: 

(1 + a)<j)" s + 2T ra J s + uol<p s = (9) 
(l-a)0: + ^ o 2 a = O (10) 

It is noteworthy that, as far as only radiation is concerned, the symmetrical mode's damping 
rate is twice the single-bubble's one, while the antisymmetrical mode is undamped. This 
feature is easily understood in terms of constructive (resp. destructive) interference between 
the acoustic waves radiated by each bubble, and parallels a well-known situation in the 
atomic physics domain (super- and sub- radiant quantum states of a couple of identical 
atoms interacting with each other through the E.M. field). From eqs. @ and (|10"D, it is 
clear that the symmetrical mode has the lower angular frequency oj s = u /\/T+~a, and the 



In our experiments A is of order lm, while d ranges from 1 to 5 cm. 
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antisymmetrical mode the higher one uo a = cjq/\/1 — a . Observe that, leaving apart the 
calculation of radiative damping, it is very easy to derive above expressions of uj s , a using 
the following trick. Let us consider the water as an uncompressible fluid (i.e. c — ► oo). The 
water displacement due to bubble z's motion simply reads: 

R 2 

u i (r,t) = d{t) f T^i (11) 

with e n = (r — r*j)/|r — r*j| = (r — r^/r*. Then, allowing for eq. (P), the overall water kinetic 
energy T is: 

= \M (el 2 + ^ 2 + 2a^'i) (12) 

where Mo = 47ri?Qpo is the effective mass of either bubble. On the other hand, the total 
potential energy V associated with the isentropic compressibility of the enclosed air reads: 

V= l -K(fi + il) (13) 

where K = ^ti^RqPq is the effective stiffness of either bubble. The Lagrange equations 
derived from L = T — V are: 

£ + a£+ufa = (14) 

<*& +£2+^06 = (15) 

which is exactly the c — > 00 limit of eqs. (0) and It is worth noticing that eqs. (]T2|) 
through dl5f) are formally equivalent to those of a system of two (L,C) electric circuits 
coupled by mutual induction with coefficient aL. In this analogy, M and K respectively 
correspond to L and 1/C, and the £j's to the electric charges of either capacitor. 
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2.3 Forced oscillation 



Let us now suppose that the above studied two-bubble system is driven by an external 
acoustic source with an angular frequency u near Minnaert's one, ojq. The phase difference 
of the driving pressures on both bubbles can therefore be neglected, since ud/c <C 1. Let 
Pei(t) be the external pressure undergone by bubble i. Motion eqs. (^) and (§) are then 
completed in: 

& + a& + r rad & + Q + wgei = -Nr ( 16 ) 



Po-Ro 

a& + £ + ru& + & + ^2 = (17) 



or equivalently: 



9F 

1 + a 

<f>l+u 2 J a = F ea (t) (19) 

with F es (t) = -(p e i(t) + Pe2{t)) I PoRo(l + a) and F ea (t) = -(p e i(t) - p e2 (t)) / p R Q (l - a). 
Solving for (p s and 4> a in above eqs. fll8|) and (]TQ|) , one gets ^i(t) and ^(t), and conse- 
quently (using eq. @) quantities pi(r,t) and P2(r,t) at any point r of the medium. At 
last, comparing external (applied) pressure p e (r, t) with the actual overall extra pressure 
p(f, t) = p e (r, t) +pi(r, t) +P2{r, t), we can experimentally measure the two-bubble system's 
response as a function of a;. In this respect (and provided that the excitation-detection 
geometry allows it), resonances are expected for uj = u s and u> — uj a . 

3. EXPERIMENTS 

Our aim is to demonstrate the existence of both above mentioned modes. From an 
experimental point of view, it turns out to be easier to implement a forced oscillation scheme 
than a free oscillation one. We therefore present the former hereafter. 



3.1 Experimental setup 



A small net (see fig.l), made up with a gauze maintained with a wire, is designed to 
catch up an air bubble in water and to fix it at any desired position without appreciably 
modifying acoustic impedance and spherical symmetry. Two such devices are used for 
studying the two-bubble system. The external driving source is a speaker and extrapressure 
p(f, t) is measured with a small microphone. A function generator, to which the speaker 
is connected, produces a c.w. sinusoidal signal with a frequency slowly swept from fi ow to 
fhigh- The signal delivered by the microphone is transmitted to a lock-in amplifier which 
compares it with the reference one (delivered by the function generator) and decomposes it 
into real and imaginary parts. Both parts can be seen on an oscilloscope and recorded with 
a computer (see fig. 2). 

In a preliminary set of experiments, without any bubble in the aquarium, the response of 
the microphone is calibrated for different speaker-microphone configurations. Two kinds of 
configurations are presented in figure 3. In figs. 3(a) and 3(b) the configuration is deliberately 
asymmetrical: the microphone is mainly susceptible to bubble 2's motion, while the speaker 
selectively drives bubble 2 (fig. 3(a)) or bubble 1 (fig. 3(b)), so that F ea {t) is nonzero: both 
modes can be excited and the associated motions detected. In fig. 3(c), the speaker is placed 
far from the bubbles; then, not only the phases, but also the amplitudes of the external 
pressures pi e {t) and p2 e (t) undergone on either bubble are appreciably equal. In such a 
symmetrical excitation configuration, F ea (t) = 0, so that the antisymmetrical mode remains 
unexcited. Observe, by the way, that since distances r x and r 2 between the bubbles and 
the microphone are equal, the latter would detect no contribution from the antisymmetrical 
mode even though it was excited (see eqs.(|3]) and (|5|): r± — r 2 and £i = — £ 2 yields pi(r 1; t) + 
p(r 2 ,t) = 0). 
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3.2 Results and discussion 



In figs. 4(a) and 4(b), the imaginary part Imp of the output signal from the lock- 
in amplifier is displayed versus the speaker frequency / for various values of the bubbles 
spacing d. Figures 4(a) and 4(b) respectively correspond to configurations 3(a) and 3(b). 
Two resonances can be made out in fig 4(a) and (though at a lesser degree) in fig 4(b). 
Observe that the sign of the signal at the higher frequency resonance is changed from 4(a) 
to 4(b), while the lower frequency one remains unchanged. This is consistent with the latter 
signal being associated with the symmetrical mode's resonance (lj s = loq/ y/1 + a < u , and 
F es unchanged from configuration 3(a) to 3(b)), and the former one with the antisymmetrical 
mode's resonance (u a = u> / V 1 — a > u , and F ea changed into — F ea from configuration 
3(a) to 3(b)). 

It is noteworthy that both resonances have appreciably the same width. This is in 
contradiction with simplified eqs.(^) and (|T0D (or ( Pf) and ([19])), in which only the radiation 
damping was considered. In fact, as mentioned in introduction, other kinds of damping 
should be taken into account: if viscous damping is absolutely negligible for such large 
bubbles, thermal damping is not (see fig.8 in fl9|D, and may be at the origin of the linewidth. 
Further discussion of this point is out of the scope of the present paper. In figure 5, we have 
plotted, for both symmetrical and antisymmetrical modes, and for Rq ~ 2mm, the inverse 
squared frequency f~ 2 (multiplied by a factor of 10 7 ) versus the inverse bubble spacing d" 1 , 
in order to get a visual check of theoretical relations: 



Although experimental points are appreciably aligned, the measured slopes are about 
40% below theoretical prediction, suggesting that coupling constant a has been overesti- 



1 1 R 1 

J? ~ 11 + 7!' d 



(20) 




1 1 Rn 1 



(21) 



Pa /0 2 /0 2 ' d 
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mated. In fact, theoretical value a = Ro/d was derived in eq.(|l2"D when integrating the 
water kinetic energy density ^Po(fj) 2 over the whole spaceQ. This inertial coupling is nat- 
urally lowered if some obstacle lies between the bubbles and consequently screens (part of) 
the water flow|. Now, this is exactly what happens in configurations 3(a) and 3(b): to be 
able to excite the antisymmetrical mode, we are compelled to insert the speaker between 
the two bubbles, thus bringing about the above screening effect. In order to check this inter- 
pretation, we performed the same experiment with configuration 3(c), and recorded, for the 
symmetrical branch of the linear fitting of fig. 5, a slope of about 90% of the theoretically 
predicted value. 

At last, it should be noted that, when blocked from the top by the net, the bubble is, 
strictly speaking, no longer spherical. As mentioned in footnote 1, such a deviation from 
the spherical shape is (almost) of no consequence, and we have used the formulas derived 
above in section 2 with R standing for the radius of the sphere of equivalent volume (i.e. 
the radius of the bubble before it is captured by the net). Nevertheless, this feature of 
our experimental setup raises the following difficulty: since, in course of motion, the fixed 
point of the bubble is no longer its centre (as implicitly assumed in the theoretical model) 
but its top, expression ([TT]) of the water displacement is no longer correct; a dipolar term 
should be added to the spherical monopolar one. As a consequence, the kinetic energy T 
derived in eq.(|12"D is modified too. More precisely, an exact calulation shows that T should 
be multiplied by a factor of 7/6 and the coupling constant a by a factor of 1 + i?g/(4(i 2 ). 
These corrections lie within our experimental accuracy. We have consequently neglected 



5 More precisely: over the whole space outside the two bubbles (the inner air's inertia being 
negligible). Nevertheless, it can be shown that the coefficient of coupling term £^£ 2 i n integral (12) 
does not depend on the bubbles radius Rq. 

6 The effective mass Mq is modified too, but at a lesser degree. 
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them. In this respect, it may be noted that the gauze in our device does not act like a rigid 
wall because the water can flow through the meshes of the net. The situation is therefore 
different from that discused in other studies |T5| , [TE| ] considering the influence of the proximity 
of a rigid boundary on the Minnaert frequency. 

4. CONCLUSION 

As a conclusion, the acoustic inertial coupling between two air bubbles in water has been 
experimentally put in evidence. Theoretical analysis shows that the two-bubble system is 
formally equivalent to a set of two magnetically coupled (L,C) electric circuits, with two 
eigenmodes, respectively symmetrical and antisymmetrical. Experimental measurements 
and theoretical predictions are in 10% accuracy agreement. 
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FIGURE CAPTIONS 



FIG.l Simple tool for capturing the bubble 

FIG. 2 Diagram of the experimental setup. In the experiment, we pump air into a tube 
immersed in water to produce the bubbles. The radii difference between these bubbles 
is small and will be neglected. (It can be shown that a small radii difference yields 
second order correction of u s and u a ) . 

FIG. 3 Different geometrical configurations 

FIG. 4 (a) Spectra of symmetrical and antisymmetrical modes in configuration 3(a). 

(b) Change of sign of Imp for the antisymmetrical mode when configuration 3(b) is 
adopted. 

FIG. 5 Linear fitting of the plot 4- ■ 10 7 {Hz- 2 ) vs. h(cm x ) for the two modes. Average of 

Js,a f ' 

the resonance frequency of the two bubbles: 1499Hz; corresponding radius: 0.217cm; 
slopes for the two fitting lines: 0.578 and —0.550 (cm ■ sec 2 ); slopes of theoretical 
prediction: ±0.966(cm • sec 2 ). 
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FIGURES 




FIG. 1 
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